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ON ONE INVERSE SPECTRAL PROBLEM RELATIVELY DOMAIN 


Y.S. GASIMOV t 


Abstract. Different practical problems, espesially, problems of hydrodynamics, elasticity the¬ 
ory, geophysics and aerodynamics can be reduced to finding of an optimal shape. The investi¬ 
gation of these problems is based on the study of depending of the functiuonals on the domain, 
their first variation and gradient. 

In the paper the inverse problem relatively domain is considered for two-dimensional Schrodinger 
operator and operator Lu = A^u and the definition of s—functiuons is introduced. The method 
is proposed for the determination of the domain by given set of functions. 

Key words: Shape optimization, inverse problems, domain variation, convex domains, support 
functions. 

Mathematics Subject Classification 2000: 31B20, jONfS, 35R30, 65N21 


1. Introduction 

One of the well studied classes of the inverse problems-is inverse spectral problems. The papers 
dedicated to the investigation of these problems traditionally focus on the construction of a 
function (potential) by given spectral data (scattering data, normalizing numbers, eigenvalues) 
and obtaining necessary and sufficient conditions providing unequivocal determination of the 
sought function. More datail review can be found in the paper [1]. 

There exists a wide class of practical problems requiring determination of the domain by 
some experimental data. For example, it is very important to hnd the domain of the plate 
under across vibrations by the quantities, which may be measured from distance [2]. There 
are different formulations of the inverse problems relatively domain for the different ceses [3-5]. 
Note that in differ from traditional inverse problems, the inverse problems relatively domain 
have some special specifications. First, these problems require to find no function, but domain. 
Second, choice of data (results of measure) sufficient for determination of the domain is also 
enough difficult problem. 

In the paper we consider formulation and investigation of one inverse problem relatively 
domain for the two dimensional Schrodinger operator that, in particular, describes the vibration 
of membrane. In the end of the work we put and solve the similar problem for the operator 
describing across vibrations of the plate. 

2. Problem setting and preliminary resnlts 

The object under investigation is the problem 

— /S.u + q{x)u = \u, X G D, (1) 

u (x) = 0, X G Sd, (2) 

where q{x) is differentiable non-negative function, satisfying condition t^q {xt) = q{x), t G 
R,0 ^ D C R^-bounded convex domain, Sd G its boundary, A-Laplace operator. 

It is known [6,p.333] that under these conditions eigenfuntions Uj{x) of the problem (^, Q 
belong to the class {D) n [D) and eigenvalues Xj are positive and may be numbered as 

< A 2 < considering their multipilicity. 
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The set of all convex bounded domains D £ we denote by M . Let 

K = £ M : Sd e , 

where is a class of the piece-wise twice continuous differentiable functions. 
Definition!. The functions 


Jj (x, D) = 


A. 


x£E^, j = 1 , 2 ,... 


are said to be s- functions of the problem o, 0 in the domain D. 
The problem is: To find a domain D £ K, such that 


( 3 ) 


Jj{x,D) = Sj{x) , X £ Sd, j = 1,2,..., (4) 

where Uj {x) and \j is an eigenfunction and eigenvalue of the problem 0-0 in the domain D 
correspondingly, Sj (x)-given continuous functions defined on R?. 

First of all we give the considerations, which led us to this formulation. 

Investigation of the dependence of the eigenvalues of the operators on the domain is an im¬ 
portant problem, as the mechanical characteristics of some systems indeed are eigenvalues of the 
corresponding operators, which may be expressed by the functionals depending on the domain 
[7,8]. One of the important steps for the investigation of the properties of these characteristics 
is calculation of the variation of these functionals relatively domain. 

But to do this we need to define the space of domains, give a scalar product and a definition 
for the domain variation in that space. 

It was shown [9] the pairs {A, B) £ M x M form a linear space with operations 

{A,B) + {C,D) = {A + C,B + D), 


{A,B) = {C,D), if A + B = C + D, 


X{A,B) = {\X\A,\X\B). 

Here A + B is taken in the sence of Minkowsky, i.e. 

A B = {tt b '. d £ A, b £ By . 
The scalar product is introduced by formula 


where 


(a, 6) = I P,{Om)dC, 
Sb 


a — {Ai,A2), b—{Bi,B2), 

Pi (x) = Pai (x) - Pa2 (x) , 

P2 (x) = Pbi {x) - Pb 2 (x) , 

Sb— unit sphere, Pd (x) = ms^{x,l) ,x £ E'^— support function of the domain D. 
The obtained space we define by ML 2 . 

For any fixed D £ M eigenvaulue Xj of the problem 0, 0 is defined as ([10], p.l82) 

Xj = inf I {u,D), {u,Up) = 0, p = l,j - I, 
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where 


f jVu (x)l^-h q (x) (x) 


I(u,D) = 


D 


dx 


f (x) dx 
D 


Thus we can consider Xj as a functional of D € K and define by Xj{D). The following formula 
is obtained (see [9], p.98) for the first variation of the functional Xj {D) in the space ML 2 


5Xj{D) = —max f \\/uj 5Pd {n {^)) ds, (5) 

Sd 

where \\/uj (a:)p = > ^(0' outside normal to Sd in the point max is taken over 

i=l ^ 

all eigenfunctions corresponding to the eigenvalue Xj in the case of its multipilicity. 

Using © the following formula may be obtained for the eigenvalues of the problem o, m 
in the domain D 


>^j{D) = ^max f \Vuj{()f PD{n{^))ds. ( 6 ) 

^ '^3 J 

Sd 

Really, let’s take Dq £ K, D (t) = t ■ Dq, t > 0. 

By Uj we define j-th eigenfunction of the problem o, © corresponding to the domain Dq. 
Then 


— ^Uj (x) + q (x) Uj (x) = XjUj (x), X G Dq. 
This relation may be written in the following equivalent from 


(f) + 


Xj (Dq) (X 
-U ' 




(I) , x€D{t). 


(7) 


Since the function 

satisfies the relation 


Uj (x) = Uj , X £ D (t) 


Auj{x) = ^Auj (^ 1 ) ( 8 ) 

from the condition Pq (tx) = q (x)and (8) one may get 

—Auj (x) + q (x) Uj (x) = ^2 Tij {x), x £ D (t). 

It shows that Auj (x) is an eigefunction, and A {t) = eigenvalue for the problem (P), 

(P|) in the domain D{t). Then using Q we can write 


(9) 


Xj {t + At) — Xj (t) = Xj (D (t + At)) — Xj {D (t)) = 

= / \'^ui^)f [PD{t+At){n{0) - PD{t){n{0)]ds + o{At), ^ £ SD(t) ■ 

Sdu) 

If support function PD(t) (x) of the domain D {t) is differentable relatively t, then dividing 
both sides of © by t we obtain 


A'(t) = -max / iVuj (x)r P^(i)(n(0)ds 


( 10 ) 


5D(t) 
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where (x) = f^Poit) (x). 
Considering this we have 


. A, {Do) 
P 


—^ max 

o Uj{x) 


Sd 


Vuo 



2 

Pd, in{0)dS, 


C(^Sd. 


Taking t = 1 from this we get 

As we see from the boundary values of the function |Vuj(x)| unquiovoccally dehne 
eigenvalue Xj. 

From ® taking into account we obtain 


J Sj iC) PD{n{^))ds = 2, j = l,2,... . (11) 

Sd 

This is the basic relation for the solving of the considered problem. 

Note. As we take s- functions as a given data, let’s consider them for some concrete cases. 
For one dimensional case 


u" + q (x) y = Xu, 


( 12 ) 


u{a) = u{b) = 0, 

where q (x) = c > 0, 0 ^ {a, b) C R, s-functions defined by ® indeed are 




Thus the expression (8) takes a form 


Jj (a), 

J, (b ). 


(13) 


Jj {h) - b- Jj (a) ■a = 2,j = 1,2,... (14) 

Let’s take a = 0, i.e. consider the problem d, d in the interval (0, b). For this case from 
® we get the following 

Consequence. All s-functions of the problem d, d satisfy to the condition 

Jjib) = ‘^, j = 1,2,... . (15) 

This formula allows to solve the inverse problem: Let the set of functions Sj (x), j = 1, 2,... is 
given. In this case the problem of finding of the domain satisfying Q is reduced to determination 
of the point b, which may be done using d- 

As noted in consequence all s-functions satisfy to the condition d, which is equivalent to 
the one condition. This condition is sufficient for finding of the point b. Really as one may get 
from d 


j,{by 

Similarly, if 6 = 0, then we have 


2 

Note that if Jj (x) = Cj, x € Sd, Cj = const, j = 1, 2,... then as follows from d they all are 
equal to each other for all j = 1, 2,... . 

In two dimensional case from d is obtained that if the functions Jj{x, D) are constant, then 
Jj{x,D) = ^^, j = 1,2,... (see [9]). 
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Now we prove the lemma that will be used later on. 

Lemmal. Let / (x) be continuous function defined on the unit shere^s. Then for any 
L>i, D2eK 

j finiO)ds = j f{n{C))ds+ j f{n{C))ds, (16) 

where Di + D 2 is taken in the sence of Minkowsky i.e. 

Di + D 2 = {x : X = xi + X 2 , Xi e Di, X 2 ^ D 2 } . 

Proof. It is known [11], that f (x) may be continuously, positive-homogeneously extended 
over all the space and presented as a limit of the difference of two convex functions 

/ (x) = lim [gn (x) - hn (x)]. (17) 

n—>-oo 

First consider 

f{x) = g{x)-h{x), (18) 

where g (x), h (x) are convex, positively-homogeneous functions. 

As is known [12] for any continuous, convex, positively-homogeneous function P (x) there 
exists the only convex bounded set D such, that P (x) is a support function of D, i.e. P (x) = 
Pd (x). The opposite statement also is true. 

It is also known that D is found as sub differential of its support function in the point x = 0 

D = dP{Q) = {I ■. P (x) > {I, x), Vx G i?"} . 

So, there exist the domains G and 77, such that 


g (x) = Pg (x) , h{x) = Ph (x) . 
Considering (|TH|). (HU) we get 

f f(n(0)ds= f Ig (n (^)) - h (n (0) ds] = 


= I PG{niO)ds- f PH(n(^))ds. 

As for anyHi, D 2 € K the following relation is valid [9] 

j Pd 2 in{0)ds= j Pdi (n(O)ds, 

from (EDI) one may obtain 


(19) 


( 20 ) 


( 21 ) 


J f{n (0) j {n (0) ds- j PD 1 +D 2 {n (0) ds. 

SD1+D2 Sh 


As PD 1 +D 2 ( x ) = Pdi (x) -7 Pd 2 (x) [12], applying (ETJ) again we get m - Lemma is proved. 


3. Main results 

Now we investigate the main problem of the work-construction of D by given set of functions 
Sj {x) , j = 1 , 2 ,.... 

Let B C he unit ball with the center at the origin and Sb- its boundary. By (pk (x), k = 
1,2,... we denote some basis in C (S'B)-space of continuous in Sb functions. These functions 
may be continuously, positive-homogeneously extended to B. It may be done as: 
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Vk{x) = 



One may check, that these functions are continuous and satisfy to the positive homogeneity 
condition 


(pk {ax) = a(pk {x), a> 0. 

Without loss of generality we can denote (fk {x) by ipk (a^)- 

Thus we obtain the set of continuous, positive-homogeneous functions defined on B. 

As we noted above each positive-homogeneous, continuous function ipj (x)may be presented 
in the form 


ipk (x) = lim (x) - (x) 

n—>oo L J 

and there exist satisfying above mentioned properties domains and H!^ such, that 


( 22 ) 


an [x) = Pg^ (x) , h^n (a;) = Ph!^ (x) ■ 

These domains we call basic domains. Substituting these into (ESI) we get 


First we consider 



(23) 


^ {x) = (x) - Phi^ (a;), (24) 

where G^and are closed, bounded convex domains. 

As n (x) € Sb, for any x G we can decompose Pd {x),x € Sb hy basic functions (pk {x) 

OO 

Pd {x) ='^ akPk {x), x e Sb, ak e R, A: = 1 , 2 ,.... (25) 

k=l 

Thus, to determine Pd (x) we have to find the coefficients fe = 1, 2,... . 

Theoreml. Let the set of functions Sj{x),j = 1, 2,... is given. Then the coefficients Ok, k = 
1,2,... of the support function of sought domain D for which is true, satisfy the equation 

OO 

^ Ak,m{j)akam = 2, j = l,2,... (26) 

k,m=l 

with coefficients 


^k,m (j) = f Sj (x) [Pcm (n (x)) - Pnm (n (x))] ds- 

Sok 

— f Sj (x) [Pg"* {n {x)) — Ph”^ (n (x))] ds. 

Sfjk 

Proof. Considering (1241) from (|25l) one may get 


(27) 


Pd{x) = '^ak{PGk{x) - PHk{x)),x G Sb- 

k=l 

The set of all indexes for which ak > 0 {ak < 0) denote by {I~)- 


( 28 ) 
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Then the relation (P|) may be written as 


Pd {x) - oikPck (x) + X] (a?) = 

k&I- k&I+ 

= c^kPck (x) - akPnk (x), x € Sb- 

k£l+ k&I- 

From last taking into account the properties of support functions [12] we obtain 


( 29 ) 


D-Y, = E - E 

kel- k&I+ kel+ k&I- 

The use of (1291) and the lemma gives 


f Sj{OPD{n{^)d( + f Sj(C)PDiniO)d^+ 


Sd 

E ( 

kGl~ 


+ 

f s,(OPD(n(Od^ = 

J Sj{^)PB{n{())d^+ 

E 

o^kSf^k E 

^kSQk 



+ 

f MO)PD(n(0)dC 


E 

kGl- 

{-Oik)Sjjk 



From this considering m we have 


f sj(OPD(n(Od^ 

Sd 


OQ 


E «fc 

k=l 


I SjiOPD(n(C))dC- 

_S^k 


- I SjiOPDiniOdC 

S[fk 


= 2 . 


Substituting here (f^ one may get (HH) with coefficients (EH). Theorem is proved. 

We assumed that the considered problem has a solution in general case. For interesting cases 
as the functions Sj (x ), j = 1,2,... are defined as experimental data, this problem always has 
a solution. The function Pb (x) is constructed by the help of the solution of dsni) , using ((2S1). 

As we noted above domain D is unequivocally defined by its support function Pb (x). Suppose 
that dsni) has the only solution providing convexity of the support function of D. 

Let’s show that the expressions j = 1,2,... for the problem (H, (0 in the con¬ 

structed by the help of this solution, using H25I) domain D indeed are s-functions. Really, if D is 
a domain in which the problem (H, m has given by formula if^ s-functions then decomposi¬ 
tion D by formulae (ESI) and making above done transformations we get the equation (EE! with 
the same coefficients. From the assumption that this equation has the only solution, it follows 
D = D. 

If (j19|) has more than one solution then the sought domain is among the ones, constructed by 
m using these solutions, providing convexity of P{x). 

This algorithm is constructed considering m- In general case when (p (x) has form (ESJ 
Ak,m (j) turns to 





Ak,m U) = lim 


f Sj (x) [Pg- (n (x)) - Ph;^ {n (x))] ds- 
S^k 

- f Sj (x) [Pg- (n (x)) - Ph^ (n (x))] ds . 
s^k 

Now let’s consider the across vibrations of the plate. 

Let P £ P^ be a domain of the plate with boundary Sd £ ■ 

It is known [2] that the function u) (xiX 2 , t) describing across vibrations of the plate satisfies 
equation 


^XlXlXlXl p ‘^^XIXIX2X2 P ^X2X2X2X2 PuJu = 0 ■ 

Assuming the process stabilized the solution - eigen-vibration is sought as 

Ld (Xl, X 2 , t) = U (Xl, X 2 ) COsXt, 

where A-is an eigen-frequency. 

Substituting this to dsni) we get 


(30) 


= Xu, x G D, (31) 

where = AA. 

For different cases different boundary conditions may be considered. The object under inves¬ 
tigation is the freezed plate with boundary conditions 

Sij 

u = 0, -^ = 0, X £ Sd- (32) 

on 

Let 


K = G M : Sd G , 

where is a class of the piece-wise twice continuous differentiable functions. 

Definition2. The functions Jj {x,D) = ^ G P^, j = 1,2,... are called s- func¬ 

tions of the problem (EU), ED) in the domain D. 

The problem is: To find a domain D G K, such that 

Jj (x,P) = Sj{x),x £ = 1,2,..., (33) 

where Uj (x) and Xj are eigen-vibration and eigen-frequency of the problem EUi-EH in the 
domain D correspondingly, Sj (x) , j = 1,2,. . .-given continuous functions defined on P?. 

In [9] the following formula is obtained for the eigen-frequency of the freezed plate under 
across vibrations 


= f \Auj{()\‘^ PDin{^))ds, (34) 

4 '^j J 

Sd 

where Pd (x) = max (l,x), x G is a support function of D, and max is taken over all 
leD 

eigen-vibrations Uj corresponding to eigen-frequency Xj in the case of its multiplicity. (As we 
see from El, the boundary values of the function \Auj (x)p unequivocally define Xj). From 
El considering El we get 
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J Sj{C)PD{n{x))ds = A, j = 1,2,... . 

Sd 

Caring out above done considerations the following theorem is proved for considered problem. 
Theorem2. Let the set of functions sj (x), j = 1, 2,... is given. Then the coefficients of the 
support function of sought domain D of the plate for which (jSl is true, satisfy the equation 

OO 

^ ^ ^k,m O') 4, j 1, 2, 

k,m=l 

with coefficients 


■^k^m {j) 


lim 

n—^oo 


/ Sj (x) [Pg- {n (x)) - Ph^ (n (x))] d,s- 
S^k 



Sj (x) [Pg^ (n (x)) - Ph^ (n (x))] ds 
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